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The magnetic flux periodicity of superconducting loops as well as flux quan- 
tization itself are a manifestation of macroscopic quantum phenomena with far 
reaching implications. They provide the key to the understanding of many funda- 
mental properties of superconductors and are the basis for most bulk and device 
applications of these materials. In superconducting rings the electrical current 
has been known to periodically respond to a magnetic flux with a periodicity 
of h/2e. Here, the ratio of Planck's constant and the elementary charge defines 
the magnetic flux quantum h/e. The well-known h/2e periodicity is viewed to 
be a hallmark for electronic pairing in superconductors and is considered evi- 
dence for the existence of Cooper pairs. Here we show that in contrast to this 
long-term belief, rings of many superconductor bear an h/e periodicity. These 
superconductors include the high-T c cuprates, Sr 2 Ru0 4 , the heavy-fermion su- 
perconductors, as well as all other unconventional superconductors with nodes in 
the energy gap functions, and s-wave superconductors with small gaps or states 
in the gap. As we show, the 50-year-old Bardeen-Cooper-Schrieffer theory 
of superconductivity implies that for multiply connected paths of such super- 
conductors the ground-state energies and consequently also the supercurrents 
are generically h/e periodic. The origin of this periodicity is a magnetic-field 
driven reconstruction of the condensate and a concomitant Doppler-shifted en- 
ergy spectrum. The robust, flux induced reconstruction of the condensate will 
be an important aspect to understand nanoscale properties of unconventional 
superconductors. 

Currents of electrons moving on multiply connected paths are modulated by an applied 
magnetic flux with a period of h/e ^\ as predicted by Aharonov and Bohm 2 . In super- 
conducting rings the order parameter responds also periodically to a magnetic flux, as Fritz 
London recognized when he analyzed the implications of a single-valued superconducting 
wave function 3 ; different condensate states, which differ by integer flux quanta, are related 
by a gauge transformation. London concluded that the flux periodicity in superconduct- 
ing rings is h/e He missed, however, a class of supercurrent carrying wave functions, 
which were identified years later 4-6 , and allowed to explain the experimentally observed 
h/2e flux quantization 7,8 . Indeed, according to the Bardeen-Cooper-Schrieffer (BCS) the- 
ory of superconductivity 9 the electronic condensate is formed by Cooper pairs, which carry 
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twice the elementary charge. However, fundamentally it is not just the pairing motivated 
substitution of e by 2e, from which the periodicity in h/2e originates, but rather the sub- 
tle requirement of the degeneracy in energy 4-6 of the two distinct classes of supercurrent 
carrying states. 

The original flux trapping experiments 7 ' 8 , which proved the h/2e flux quantization in 
superconductors, as well as the later experiments 10-12 were considered a manifestation of the 
formation of Cooper pairs in the then known conventional superconductors. The discovery 
that magnetic flux changes the magnetization of YBa 2 Cu 3 7 _ ( 5 rings with a periodicity of 
h/2e was similarly argued to provide the evidence for Cooper pairs also in high-temperature 
superconductors 13 . 

Does, vice versa, the existence of Cooper pairs or the h/2e flux quantization necessarily 
imply an h/2e periodicity of the energy or the current in superconducting loops? In fact, the 
h/2e periodicity requires that multiply connected superconductors threaded by a flux nh/2e 
are degenerate in energy for different integers n. In superconducting s-wave rings or hollow 
cylinders with inner diameter d this degeneracy occurs if d ^> £, where £ is the coherence 
length 4-6 . In the opposite regime d < £ the discrete quantum nature of the electronic states 
in the ring matters and the energies at half-integer and integer flux quanta are generally 
different; correspondingly the superconducting behavior is only /i/e-periodic (see Fig. 8-8 
in Ref. 14). This behavior should be observable, possibly in Al rings with d < £ = 1.6 fxm. 

The oscillation period of energy and currents in superconducting rings is therefore not 
always h/2e. In fact, as we report here, the BCS-theory strictly predicts that for rings of 
superconductors with nodes in their gap functions, such as the high-T c cuprates, Sr 2 Ru04, 
or the heavy-fermion superconductors, the ground-state energy is generically h/e periodic. 
For all these superconductors, the states that yield the BCS-condensate state also include 
current-carrying states with energies close to the Fermi energy E F . As a result of the 
magnetic-field driven change of occupation of these states and the concomitant reconstruc- 
tion of the condensate, the superconducting rings develop an h/e periodicity of the super- 
current. 

The flux periodicity in mesoscopic loops of rf-wave superconductors is contained in the 
solution of the Bogoliubov-de Gennes (BdG) equations 15 for the pairing Hamiltonian: 

(ij) ,tr (ij) 
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Figure 1 | Real-space representation of a square loop with a typical electronic prob- 
ability density |\I/| 2 of a single state in the condensate. The figure displays an eigenstate 
of the d-wave pairing Hamiltonian, calculated for a square-loop with 80x80 lattice sites with a 
pairing interaction of 0.3t. The hole in the center has a size of 28x28 unit cells. To enhance the 
contrast of the complex pattern, the special color code shown on the right is used and the discrete 
lattice points are smoothly interpolated. 

The operators Cj a (cJ CT ) annihilate (create) an electron on lattice site j with spin a =f, |; 
t is the hopping matrix element between nearest neighbor sites, (fij = 2ne/h J. A(r) ■ dr 
is the Peierls phase factor, and A is the vector potential of the magnetic field. The order 
parameter of the superconducting state A^, denned on the links between neighboring sites 
with phase factors appropriate for <i-wave symmetry. 

Figure 1 displays the probability density of the wave function for a state with energy close 
to Ep on a square loop, whose edges are oriented parallel to the [100] and [010]-directions, 
respectively. The (i-wave loop eigenstates are obviously far more complex than the angular 
momentum eigenstates of a one-dimensional circular ring (cf. Ref. 14), and the current flow 
in this loop can only be evaluated numerically. Nevertheless, also a qualitative discussion 
allows insight into the underlying physics. 

To assess the global quantities, viz. energy and current, the evolution of the eigenenergies 
with magnetic flux has to be calculated. The eigenstates with energies below Ep form the 
ground-state condensate (Fig. 2). Only flux values $ between and h/2e are discussed, 
because all quantities are either symmetric or antisymmetric with respect to flux reversal 
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Figure 2 | Energy spectrum of the cZ-wave BCS model. The eigenenergies in the gap region 
are shown for a square 40x40 loop with a hole of 14x14 unit cells and pair interaction 0.3i as a 
function of flux <]? (in units of h/e). The energies are given in units of the superconducting order 
parameter Ao at $ = (Ao ~ 0.22t). The superconducting condensate consists of the states below 
Ep = (red lines). Reconstruction of the condensate takes place near $ = ±(2n + l)h/Ae, where 
the eigenenergies jump abruptly. 

$ — > — $. Two clearly distinct regimes are found: the flux intervals between and h/Ae and 
from h/Ae to h/2e. 

Up to $ = h/Ae the supercurrent J generates a magnetic field which tends to reduce the 
applied field. This is achieved by a continuous shift of the eigenenergies in the condensate. 
At $ = 0, pairs of states with opposite circulation compensate their respective currents, 
thus J = 0. The well separated states at $ = in Fig. 2 are the states in the vicinity of the 
nodes of the mesoscopic d-w&ve superconductor. At energies further away from Ep, the state 
density is higher; these are the states near the maximum energy gap that provide most of 
the condensation energy. For $ > 0, the energy of the states with orbital magnetic moment 
anti-parallel (parallel) to the magnetic field is increased (decreased). Correspondingly the 
supercurrent, which is carried by these states, depends on the details of level crossings and 
avoidings. The main contribution to the supercurrent arises from the occupied levels closest 
to Ep, because the contributions from the lower-lying states tend to cancel in adjacent pairs. 

As the highest occupied state shifts with increasing flux to lower energies, the current 
in the square loop first increases for small $ (Fig. 3), then decreases, when the highest 
occupied level with an orbital moment opposite to the applied magnetic field starts to 
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Figure 3| Flux dependence of energy and current. Total energy (E($>) — E(0))/E(0) (a) 
and total circulating current J (b) for a square 40x40 loop with a hole of 14x14 unit cells and 
pair interaction 0.3 1 as a function of flux in units of h/e. J is given in units of et/h = 6 x 10~ 5 A 
for a typical choice of t = 250 meV. There is a clear difference between condensate states with an 
even and an odd winding number q of the order parameter, reflected e.g. in the deformation of the 
q = 0-parabola. The overall ^-periodicity for E and J is h/e. 



dominate. With increasing flux this state approaches Ep. For s-wave rings this "Doppler 
shift energy" (cf. Ref. 15) corresponds to the critical value of the superfluid velocity, for 
which the indirect energy gap closes. For <i-wave loops, the order parameter is protected by 
the numerous states that form the "lobes" of the (i-wave gap parameter. 

For d-wave loops and rings with other unconventional order parameter symmetries, the 
states in the vicinity of the nodes evolve with increasing flux as in small gap s-wave rings. 
They do not necessarily cross Ep (Fig. 2) due the hybridization of the respective states 
above and below Ep. Nevertheless, a state with one direction of current is replaced by 
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a state of opposite direction (Fig. 4). The current carrying states of the condensate are 
thereby continuously changing near the extrapolated crossing points. As a consequence the 
energy "parabola" centered at zero flux is different from the ground-state energy parabola 
centered at $ = h/2e (Fig. 3). The deviation from a parabolic shape near zero flux is 
due to the evolution of the near-nodal states; the vertical offset of the energy minima at 
$ = nh/e results mostly from the flux dependence of the states near the maximum value of 
the anisotropic gap. 

For a flux value near h/Ae the condensate reconstructs. The superconducting state beyond 
h/Ae belongs to the class of wavefunctions introduced by Byers and Yang 4 . Remarkably, in 
the flux interval from h/Ae to h/2e, a full energy gap exists also for rf-wave superconductors 
(Fig. 2). Here the circulating current enhances the magnetic field; the paramagnetic moment 
of the current is parallel to the field. The resulting energy gain is responsible for the 
field-induced energy gap. This reconstruction of the condensate is the origin of the h/e 
periodicity in energy and current. Intriguingly, for superconductors with unconventional 
order parameter symmetries also larger loops (d ^> £) are h/e periodic. 

The numerical solution of the BdG equations with a self-consistency condition for the 
order parameter is adequate for ~ 15 nm rings. However, to examine systems of micrometer 
size, the nodal states have to be described using a continuous gapless density of states. The 
flux induces a Doppler shift which modifies the states and alters their occupation near Ep, 
thereby causing an h/e component of the current J. While the h/2e component of J oc 1/d, 
the h/e component decreases with 1/d 2 (see Appendix B). In quantitative agreement the h/e 
component which, as compared to a ring of the size shown in Fig. 1, reduces by a factor of 
60 for a corresponding ring of 1 /xm size, measured by the weight of its Fourier peak. Using 
typical parameter values for a YBCO ring of 1 jum size, the ratio of the h/e versus the h/2e 
component remains in the percent range. The frame width w of the ring has little influence 
on the weight of the h/e component for the loops with w smaller than the penetration depth 
A. A similar behavior is also shown by loops with w > A, because only states that result in 
the current-transport channels within A affect significantly the h/e component. 

Our calculations show that while changes in geometry, the number of transverse chan- 
nels and elastic scattering by impurities modify the J(<&) characteristics in detail, they do 
not eliminate the h/e component. As long as the single particle states are well defined, 
also electronic correlation effects, which are responsible for the renormalization of states 
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Figure 4| Current distribution in a square loop of 40x40 lattice sites. The current 
expectation value of the occupied state closest to Ep = is shown for flux $ = 0.17 h/e (top 
panel, the state is marked with 'a' in Fig. 2) and for $ = 0.21 h/e (bottom panel, marked with 
'b' in Fig. 2). The color encodes the projection of the current onto a square path around the loop 
whereby red presents a counterclockwise and blue a clockwise circulation. The maximal current is 
•/max = 0.15 et/H for a and 0.13 et/h for b. The current distribution of each of the two states has 
strong spatial variations and does not fulfill the continuity condition which, however, is restored 
for the total current. 

and of coupling parameters, are not expected to bear a strong influence on the discussed 
phenomena. 

The robust, magnetic-flux induced presence of currents that flow opposite to the main 
screening currents affect many properties of unconventional superconductor. Of particular 
importance are a resulting enhancement of the London penetration depth and a weakening 



8 



of the rf-shielding. The h/e periodicity of the supercurrent is a fundamental property of 
loops formed by unconventional superconductors. 
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APPENDIX A: NUMERICAL METHOD 



To investigate ring geometries with finite width, we self-consistently solve the Bogoliubov 
de Gennes (BdG) equations on the square frame shown in Fig. 5 for the Hamiltonian 



(ij) 



where 



V 



ho 



(Al) 



(A2) 



is the order parameter defined on the two neighboring lattice sites % and j. The pairing 
interaction strength is V and appropriate phases for rf-wave pairing: A iji+£ = — Aj i+ ^ are 
implicitly incorporated. A magnetic flux is represented by the Peierls phase factor ipij = 



2ne fJ 
h 



p A(r) ■ dr. We choose a vector potential of the form A = <p(y, —x)/(2irr 2 ), yielding 
a flux threading the hole with no magnetic field penetrating the superconductor, where 
(f) — $ e/h measures the flux in units of h/e. The Hamilton operator (Al) is diagonalized 
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Figure 5 | A two dimensional discrete lattice for a square frame with open boundary conditions. 
The figure shows the standard geometry for a system size of 40x40 lattice sites with a centered 
hole of 14x14 lattice sites, on which calculations were typically performed. 



by the Bogoliubov transformation 



EE 
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^2 [«niTn 



(A3) 
(A4) 
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where rf ns and 7„ s are creation and annihilation operators for fermionic Bogoliubov quasi- 
particles. The coefficients u n and v n have to fulfill the equation 

(A5) 

where the operators Hq and A act on the "single particle" wave functions u n and v n as 

H u ni = -ty^e'^M,^ - fJ,U ni , (A6) 
j 

Av ni = ^AijVnj, (A7) 
3 

and denotes the sum over all nearest neighbor sites of i. The order parameter is 
calculated self-consistently from 



Ai > = ^ E K«< y + u nj v* m ] tanh 



where the sum runs over the positive eigenvalues E n only and T is the temperature. The 
current density from lattice site i to j is 

J ^ = E ('V>' ; ''- - c\ jCia e*») (A9) 

s 

= -4| £ Im [( Mn ,<J(E„) + v* nj v ni (l - f(E n ))) e**] ; (AlO) 

n 

/(E) = l/(l + e s / fesT ) is the Fermi distribution function. The self-consistent solutions of the 
BdG equations on the square loop are characterized by the winding number q of the phase of 
the order paraneter Ay around the loop. For a fixed value of flux 0, ground-state solutions 
in different g-sectors are found by choosing suitable starting values for the iterations in the 
self consistency loop. 



APPENDIX B: MULTI- CHANNEL MODEL FOR LARGE D-WAVE RINGS 

Since the numerical method outlined above is not suited for calculations on loops of larger 
size say in the jum range, we use a multi-channel ring model which allows for an analytic 
calculation. A superconducting ring is thereby composed from many one dimensional (ID) 
loops with different radii. Each loop represents one current channel; the properties of the ring 
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are obtained by integrating over its thickness. In ID, the only spin-singlett pairing symmetry 
possible is s-wave pairing. To obtain the current characteristics of a rf-wave loop, we use the 
following sceme: We obtain the supercurrent in the loop through an energy integration over 
the current contribution of all occupied eigenstates of the BCS Hamiltonian. In a circular 
loop (with no hybridization), the Doppler shift of the eigerenergies is a linear function of the 
flux independent of the pairing symmetry. The only way in which the symmetry influences 
the supercurrent is through its characteristic density of states (DOS). We therefore perform 
the following calculations for a ID s-wave loop and obtain the d-wave supercurrent for a 
quasi ID channel by inserting the (Doppler shifted) d-wave DOS (Fig. 6) into the final 
energy integration for the supercurrent (Eq. (B15)). 



1. Superconducting State in a ID s-wave loop 

We describe the kinetic energy of the electrons on an individual flux threaded ring with 
N discrete lattice sites 

^o = 5^ e k-4>cis c ks, (Bl) 

ks 

by the tight binding dispersion 

e k ^ = -2tcos f k -^- ) . (B2) 



R 

ek-<f> is the energy of a single particle state with angular momentum hk with fceZ. The 
radius of the ring measured in units of the lattice constant a is R = N/2n. The BCS pairing 
Hamiltonian has the form 

H = H ° + E [A£(g)c_ fc+(?i c feT + A fc (g)4 T cL fc+ J , (B3) 

k,q 

where 

Ak(g) = Y, V f~ K c fc'i c -fc'+<?t) ~ ( c k>i c ~k>+ q i)} (B4) 

k 

is the superconducting order parameter and q G Z its winding number. For a perfectly 
circular ring geometry, the winding number can be identified with the angular momentum 
hq of a Cooper pair. Choosing the pairing energy Vkw = V independent of k and k! leads 
to pairing in the s-wave channel, which is the only possibility in one space dimension. 
Since we are interested in low temperature properties, we assume that the superconducting 
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condensate in its ground state is characterized by the quantum number q((j>), which changes 
its value as a function of flux whenever the total energies for two different q- values become 
degenerate. Up to finite size effects, the g-number of the ground state changes to the next 
integer whenever crosses the flux values (2n — l)/4, n e Z: 

,(0) = ^+^), (B5) 

where for positive (negative) x, hit (a;) is the largest (smallest) integer number equal or 
smaller (larger) than x. We therefore choose an ansatz for A k (q) of the form 

A fc (g') = 5 q(M A k , (B6) 

with arbitrary /^-dependence of For s-wave pairing Afc = A is constant and Ak(q) = 
A(q). With this ansatz the diagonalization of the Hamiltonian (B3) leads to the energy 
spectrum 

E±(k, 0) = £fc - ~ 2 £ ' fc ^ +g ± ^A 2 + (^z± ± e ^+^ \ (B7) 

The energies E±(k,(f>) shift with flux and for ^ n the particle-hole symmetry of the 
spectrum is broken. Near the Fermi energy, in different g-sectors, the Doppler shift e(0) of 
the eigenenergies is found by expanding E±(k, 0) in 0, leading to e(0) = ±(0 — q/2) 2t/R + 
0((4>/R) 2 ). If therefore A(q) < e(0 = 1/4) = t/2R, where the condensate changes q form 
to 1, the indirect energy gap closes and the occupation of states changes. For > 1/4, 
the pairing of electrons in states with total angular momentum q ^ according to Eq. (B4) 
becomes favorable. In one dimension, the gap closes exactly at the depairing velocity of the 
condensate beyond which no self-consistent solution of the order parameter exists 16 . 

2. Current in a d-wave loop 

In this section, we first derive a expression in form of an energy integration for the 
supercurrent in a s-wave loop which is then transformed into a rf-wave loop as described 
above. 

In the nearest-neighbor tight binding model for a ID ring, the current is given by 

J(R) = IY1 MR)** = ^ E ^f^), (B8) 

ks ks 
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where J k = de k -$jdk is the group velocity of the state with angular momentum hk and the 
occupation probability n ks is obtained using a Bogoliubov transformation: 

n ks = (clc ks ) = £ f ( , e ^ + £ -^ ±^ = + a ) f(E a (k,<p)). (B9) 

Eqs. (B8) and (B9) are a closed-form solution for the total current in a superconducting flux 
threaded ring. The sum over k has to be computed numerically, though. As shown below, 
the expansion of Jk{R) in powers of <f)/R, provides a <fi independent contribution which is 
paramagnetic for q = and diamagnetic for q = 1 plus a contribution linear in 0, which is 
diamagnetic for q = (Meissner effect) and paramagnetic for q = l 17 . 
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Figure 6 | Scheme for the density of states of a ti-wave superconductor for a finite flux |0| < 1/4 
{q = 0) (a) and for 1/4 < <f> < 3/4 (q = 1) (b). The energies are Doppler shifted to higher (red) 
or lower energies (blue). This results in a double-peak structure of the coherence peaks and for 
q = in an overlap of the upper and lower band in the region — e(<p) < E < e((f>) ( 18 \ States in 
the upper band become partially occupied. For q = 1 there is a gap of width 25i = 2o"i(0). The 
black line in b represents the density of states for eft = 1/2. 




14 



In the following we restrict the discussion to the flux interval — 1/4 < < 1/4 where 
q = in the ground state. We assume that R >> 1 and expand e ±fc _^ and J k in <p/R for 
k > as: 

2* 



e± — e±fc-(A — e fc ± 



J ± (i?) = J ±k (R) 



2te 
R~ 



1 



) +O((0//2) 2 ), 
e fc 



R2t 



(BIO) 
(Bll) 



To leading order, the quasiparticle energies in the superconducting state become 

E ± = E ± (±k, 4)=E± |^1 -(|) 2 + 0{{<P/Rfl (B12) 

with £ = sje\ + A 2 if e fc > and £ = -a/^! + A 2 if e k < 0. In the vicinity of the Fermi 
energy Ep, this simplifies to e± = efc ± e(0) and .E± = 22 ± e(0). Converting the sum over k 
in Eq. (B8) to an integral over the normal state energy e k , the total current becomes 

J(R) = J AeD{e)n + {e)J + {R,e) + J de D(e) n_(e) J-(R, e) + 0{{4>/R) 2 ), (B13) 



where 



a=±l x v 7 



(B14) 



and D{e) is the DOS of the normal state. With e = ±^ E 2 — A 2 we rewrite Eq. (B13) as 
J(R) 



f2t 

N I dE 

-it 



J2 l -(a + D s (E))f(E + 5(<P))J-(E) 
«=±i 2 

+ E J(«+^))/(^-w)j + (s) 



Q = ±l 



(B15) 



where we assume D(e) = Nq constant in the vicinity of Ep and D S (E) is the DOS in the 
s-wave superconductor: D S (E) = \E\/^E 2 - A 2 if \E\ > A and D S (E) = if \E\ < A. At 
T = the current can be separated into two contributions J(R) = Ji(R) + ^(-R), where 
Ji contains all the contributions from the states which are below E F within the interval 
— 1/4 < < 1/4 representing the standard supercurrent and J 2 contains the additional 
contributions from J(R) which appear if e((f>) > A. One finds: 



MR) = -N Q 



2e 

hR 2 ' 



-e(4>) 



dE\E\ = -No 



2t 



h 



MR) = No— / dED s (E) + 
tin j_ e (0) 



3 
R* 



(B16) 



(B17) 
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Figure 7 | The current in a thin d-wave ring as a function of flux $ (in arbitrary units). Shown 
is the result of the multi-channel model for the characteristic value 2i/(Ao-R) = 0.4. For —h/4e < 
$ < h/4e, where q = 0, the current is reduced by a contribution proportional to <£ 2 , whereas it is 
strictly linear in <J> otherwise. This gives rise to an overall current periodicity of h/e. 



where the upper integration boundary in Eq. (B16) is extended to zero. 

We replace now D S (E) in Eq. (26) by the DOS D^(E) of a d-wave superconductor as 
shown in Fig. 6 a. For finite flux 0, all energy levels are shifted according to the magnetic 
moment of their current; this results in a Doppler shift of the coherence peaks 18 . In the 
relevant regime Ao > e (</>), it is sufficient to approximate Dd(E) ~ \E\/Aq (Fig. 6 a) and 

N 2te f eW 



A hR 



dE E 



N e 2 (2t 



e(4>) 



Anh 



R 



(B18) 



The total current J(R) = Ji(R) + J 2 (R) of this channel becomes 



J{R) = -N 



h 







6 2t 1 


(I) 2 


1 - 





(B19) 



The normal state DOS, iVo = Rj (2t), is itself a function of R. The total current J for q = 
in a ring of finite thickness D and inner radius i?< is obtained from 



J 



R< 



dR(Jx(R) + J 2 (R)). 
d 



1 



2t D + 2i? < 



' N °h^ 2t ^ R < (R < + D) f r A 2R < {R < +D)_' 



(B20) 
(B21) 



In the limit of thin rings (i?< ^> D), we introduce d = 2R K in units of the lattice constant 
a and find that the ratio 



Ji Ao d 

shows the same power law in 1/d as for a single channel. 



(B22) 
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For q — ±1, an energy gap A(q) > e(0) persists for all 0, thus a flux induced effective gap 
25i(4>) is present, as shown in Fig. 6 b. Therefore calculations as above are valid also for this 
flux window, however J^{R) = 0. Only the standard supercurrent component contributes 
to the total current with J(R) = J\{R) for each channel (Fig. 7). 

Because J2 is finite for even q but zero for odd q, whereas J\ is identical for all q, we find 
that Ji is periodic with h/2e and J2 with h/e. The result in Eq. (B22) implies that the ratio 
of the h/e and the h/2e Fourier component of the total current scales with the inverse ring 
diameter. 
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